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I. THE MODEL OF THE BACTERIAL FLAGELLAR MOTOR 



In this Supporting information we provide background information on our model of the bacterial flagellar motor. 
We also derive the analytical solution of our coarse-grained model of the switching dynamics and explain the hybrid 
stochastic algorithm used for the simulations. 

u 

d 

. A. Stator-Rotor interaction 

cr: 

The model for the stator-rotor interaction is discussed in the sections The stator-rotor interaction and The rotor 
switching dynamics of the main text. The model is based on the model of Oster and Blair and coworkers [H, H], but 
^ ■ extended to include the conformational transitions of the rotor protein complex. Here, we discuss aspects of the model 
' that are not discussed in the main text. But, for completeness, we also give the main equations already presented in 
00 . the main text. 

In our model, each stator-rotor interaction is described by 4 energy surfaces, C/J. , with the subscript Sj = 0,1 
, denoting the conformationals state of stator protein j and the supersript r = 0, 1 denoting the conformational state 
of the rotor (clockwise or counterclockwise). We assume that the stator proteins are fixed by the peptidoglycan layer 
' and that only the rotor complex moves. The equation-of-motion of the rotor is then given by 

oo ; 

9: ^duuo,) 



in^--T.^0^ + P^ + m.it). (1) 



' Here, 7r is the friction coefficient of the rotor; Ug.{6j) is the free-energy surface shown in Fig. 2 of the main text. 



where 9j ~ 6^ — 6g^, with 0r the rotor rotation angle and the fixed angle of stator protein j; r]^{t) is a Gaussian 
white noise term of magnitude ^/2kBT7H ', Ns is the number of stator proteins. The torque Fj^ denotes the external 
load. As discussed in [ll, 0, 0], for the system studied here, the torque-speed curves under conservative load and 
viscous load are identical. However, as discussed in the main text, the type of load does markedly affect the CW ^ 
CCW switching dynamics. 

The transition (or hopping) rate for a stator protein to go from one energy surface to another depends upon the 
free-energy barrier separating the two surfaces. We make the natural phenomenological assumption that the hopping 
rate depends exponentially on the free-energy difference, in a manner that obeys detailed balance. Furthermore, 
following Blair and Oster and coworkers, we assume that the access of the periplasmic protons to the stator-binding 
sites is triggered by a rotor-stator interaction [H, 0|- This yields the following expression for the hopping rates: 

K.^sr (Oj) = fcou.(0,) exp[A[/,,, {0j)/2], s, s' = 0, 1. (2) 

Here, /cq sets the basic time scale, and AUss'{&j) ~ Us'{9j) — Us{6j). The function w{9j) describes the proton hopping 
windows (see Fig. 2 of the main text), which refiect the idea that the ion channel through the stator is gated by the 
motion of the rotor. 
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FIG. 1: The torque-speed relation (symbols)as predicted by the model used here, which is based on the model of Oster et 
al. The blue symbols correspond to the regime in which the load pulls the motor in the backward direction; the red 

symbols correspond to the regime in which the load pulls the motor in the forward direction, which is the scenario right after a 
switching event. The line shows that, to a good approximation, the speed a; as a function of the motor torque tm, is given by 
the speed as a function of the conservative load tl, according lo uj — {k^^^ — k^^^)'K /26, where the hopping rates in the forward 

(+) and backward (-) are given by k^^^ ~ Jq^^'^^ d6Ps{9)k^{9), and the stationary distribution Ps{S) is approximated by the 
equilibrium distribution Pcq{0) oc fi \U s{0) +tlS]) (See Ref. 0]). The parameters used in the simulations are shown in 

Table U 



The rotor complex is modeled as an MWC model [5j, which means that all the rotor proteins switch conformation 
in concert. This leads to the following expression for the instantaneous switching rate: 

¥-^'{{6,]) = fco exp[AC/'-'-'({0,})/2], r, r' ^ 0, 1, (3) 
where At/'"''' = Y^^li ^I'ji^j) ~ ^Iji^j)- The average, effective switching rate is given by 

fc^^^t'ch- j de^p{e^)K-;^' [{9,}), (4) 

where P{Oji) is the stationary distribution of the rotor's position. The instantaneous switching rate k^T^^ ({^j}) 
does not depend upon the load. Indeed, in our model, the load does not directly affect the probability that the 
rotor proteins switch conformation. In this respect, the mechanism that we propose differs fundamentally from that 
often used to explain the force dependence of processes such as protein unfolding and molecular dissociation in 
that mechanism one assumes that the reaction coordinate can described by a single order parameter, and that the 
force directly couples to that coordinate, changing the relative stability of the two (meta)stable states, as well as the 
location and stability of the transition state separating them. In our model, the propensity for the rotor to switch 
depends on interactions with the stator proteins. Consequently, the reaction coordinate for switching depends not 
only on the coordinate describing the conformational state of the rotor protein complex, but also on the coordinates 
describing the positions and the conformational of the stator proteins. While the load may change the free-energy 
landscape in the direction describing the conformational state of the rotor, we assume that the load only couples to 
the rotation direction of the rotor. The load thus changes the steady-state distribution of the rotor's position relative 
to that of the stator proteins, which in turn affects how often during their motor cycle the stator proteins favor one 
conformational state of the rotor protein complex over the other. In other words, while increasing the load does not 
change the instantaneous switching rate k^~*^ it does shift P{0b.) to positions 9j{ where kl~*^ ({^j}) is large. 

This is the principal mechanism that, according to our model, makes the effective switching rate ^g^t'ch sensitive to 
load and speed. 

The load In the experiments of Korobkova et al. the motion of the flagellum is visualized via a latex bead connected 
to the flagellar filament The bead exerts a force on the rotor protein, which, effectively, tilts the energy surfaces 

shown in Fig. 2 of the main text. When a) the connection between the load and the motor is soft, b) the dynamics 
of the motor is much faster than that of the load, and c) chemical transitions lead on average to a fixed translation 
distance of the rotor, as in the current model, then the torque-speed curves under conservative load and viscous load 
are identical P, H, [I]- However, as discussed in the main text, the type of load does markedly affect the CW CCW 
switching dynamics. 

The stator proteins Resurrection experiments suggest that in vivo the number of stator proteins is around 8 — 12 
fll| . At high load, the stator proteins act cooperatively, and the motor speed increases with the number of stator 
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proteins ; the model of Oster and coworkers describes this observation [![ . Recent experiments by Yuan Berg show 
that near zero external load, the speed is independent of the number of stator proteins [l^l- The model of Oster and 
coworkers can reproduce this behavior if the stator proteins are connected to the rigid framework of the cell wall via 
very soft springs. However, to generate a speed that is independent of the number of stator proteins, the springs have 
to be made so soft that they stretch a distance of order 10 nm, which, as Yuan and Berg point out, seems unlikely p^ . 
We therefore focus here on a motor that has only one stator protein that is rigidly connected to the cell membrane. 
This motor has a lower maximum torque than a "wild- type" motor with 8-12 stator proteins, but this is not critical, 
since we take a rather small bead (see table U]); in essence, to a good approximation, all the torques in our model could 
be scaled by the number of stator proteins. More importantly, our model correctly predicts the maximum speed of 
about 300 Hz, as recently observed by Yuan and Berg [l3|, and its torque-speed relation exhibits the distinc knee at 
a speed of about 250 Hz (see Fig. [T|). This model thus captures the effect of the dynamical interplay between the 
torque and the speed, and on the other hand the switching dynamics. The maximum speed is particularly important, 
since that, together with the total change in the winding angle of the flagellar filament upon a motor reversal, directly 
affects the characteristic switching time. In future work, we will investigate the effect of the number of stator proteins 
on the switching dynamics. 

The parameters of the rotor-stator model, which were mostly taken from [l| , are summarized in table IH 



B. Elasticity of the flagellar filament 



We assume that the free energy of a flagellar filament in a given polymorphic state m is quadratic in the curvature 
K and torsion r: 

U^,{t, k)/L = ]^EI{k - Kraf + ^tiJ{T - r„0^ (5) 

where L is the contour length, E and fi are the Young's and shear moduli, / and J are cross-sectional moments, and 
Km and Tm are, respectively, the spontaneous curvature and torsion of the filament in state m. The curvature and 
torsion are functions of the height z and the winding angle 6: 



Oz 



K{d,z) = , (6) 



rie.z) = j^. (7) 



We assume that at each instant the length of the filament has relaxed to its steady state value Zcq{9), obtained as 

dz 



a solution of the equation ^^("^(^-f )'''(^'^)) = This allows us to climate z and express the "torsional" energy as a 



Parameter 


Value 


Description 


d 


2tv rad/26 


Potential periodicity 


pmf 


152 mV 


Proton-motive force 


AG 


11.8 ksT 


AG = 2e X pmf 


ko 


3.5 10* s"^ 


Ifopping prefactor 


di 


0.05 d 


Position potential maximum 


d2 


0.1 d 


Position start power stroke 


da 


0.9 d 


Center of hopping window 


d4 


0.2 d 


Width of hopping window 


hi 


25 ksT 


Height potential maximum 


h2 


10 ksT 


Height start power stroke 


Fm = -h2/{d - 


-d2) 46A:Brrad"^ 


Force motor during power stroke 


7m 


1.7 10"-' s rad"2 


Friction coefficient motor 


ko 


0.3 s"^ 


Switching prefactor 


7L 


0.51fcBr s rad-2 


Friction coefficient load 



TABLE I: Parameters for the rotor-stator model as used in the simulations (see also 
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Parameter Value Description 

dm — dm-i 150d Spacing of wells 

ke IksT/rad^ Stiffness 

N 10 Number of wells, 

fco lO^^s^^ Jumping prefactor 



TABLE II: Parameters describing the flagellar filament. 



function of the winding angle 0: 



ul[e) = ul{T{e,z,^{o)),n{e,z,^{e))), 



(8) 



The function Uf^iO) is, in general, a complicated fonction of 6] nevertheless in the limit of equal bending and twisting 
stifnesses {EI = fiJ) the torsion potential corresponds to a linear elastic potential 



Uli9) = -keie - 9^)^ 



(9) 



where kg 



^ and 9„ 



The eperimental data of Darnton and Berg [l3| confirm that the approximation 



EI ~ fiJ is valid and therefore, locally, the potential energy guiding the dynamics of the twisting angle 9 has a 
simple linear elasticity form with elastic constant — 100 pN nm/rad^ (obtained from EI = fj,J = 3.5 pN fj,m^ and 
L — 7.6, 19.6 ^m as in [l3|)- As described in the main text, we assume that the potentiall wells are equally spaced, 
are of the same depth and have the same curvature. Clearly, these assumptions could be relaxed by allowing, e.g.^ 
the normal state to be more stable and to have a higher stiffness. 

Motivated by the observations of Darnton and Berg [l^ , we assume that the transition from one polymorphic state 
to another is an activated process, with a rate constant given by 



fc„.(0) = hexpiiU^M - Ul,{9))/2] 



The equation-of-motion for the bead is given by 



(10) 



[ID 



Here, 7l is the friction coefficient of the bead, and 77L is a Gaussian white noise term of magnitude y/2kBT^h- 
The parameters of the model are given in table |TT1 



II. COARSE GRAINED MODEL OF THE SWITCHING DYNAMICS 



We model the switching dynamics as a memoryless two-state system with switching-time distributions tp+ (t) switch- 
ing from CW to CCW and V'-(t) for switching from CCW to CW: 

CW ^ CCW (12) 

i'-it) 

Lack of memory means in this context that the probability to switch from one state depends only on the time since 
the transition to that state - the system forgets everything before the last transition. 

The switching-time distribution is related to the switching rate or switching propensity (the switching probability 
per unit amount of time) ka{t) as 

7^„(t) = fc„(0e-^o'° (13) 

One important characteristic of the stochastic trajectory of the system is the correlation function C(<) of the 
characteristic function x{t)- 



C{t) = (x(t)x(o)) - {xf- 



(14) 
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FIG. 2: A) A piecewise linear model of the switching-propensity function k{t). B) Computed power spectrum S{ui). 



We take x(0 = 1 if the system is in the CW state and xi^) = otherwise. From the ensemble of all possible 
trajectories only the ones that are in the CW state both at time zero and at time t contribute to the correlation 
function at time t. Therefore, one can write the correlation function as 

C{t) = [P(CW, t; CW, 0) - F(CW, oo; CW, 0)] P(CW, oo; CW, 0), (15) 

where P(CW, t; CW, 0) is the probability that a trajectory is in the CW state at time t given that it starts in that state 
at time zero. Using a well established result in the theory of (alternating) two-state, memoryless renewal processes 
(see [l^ . Chapter 7) one can express this quantity in the Laplace domain as: 

Pi^) ^Ul-^). (16) 
^ ' z\ ztcwj ^ ' 

Here, P{z) is the Laplace transform of P(CW, t; CW, 0), 

/•oo 

P{z)^ P(CW,t;CW,0)e-"*dt, (17) 

"'0 

G{z) is a function that depends on the Laplace transformed switching-time distributions, 

Giz) . i^-Mm^^U^))^ (18) 

(i-V-(^)V'+(^)) 

and tew is the average residence time in the CW state. The probability to be in the CW state is given by the average 
residence times as 

P(CW,oo;CW,0) = -^5^^ . (19) 

lew + tccw 



Also, using the properties of the Laplace transform one has 



C{z) = 

tew + tccw 



P{z) 



z{tcw + tccw) 



(20) 



Once we have the correlation function, wc can compute the power spectrum using the formula 

S{uj)^2 C{t) cos ujt = C{iuj) + C{-iuj), (21) 
Jo 

such that 

S{u;) = — ^— [Giicu) + Gi-iLu)] . (22) 

^-[tcw + tccw) 

In general, an analytical formula for the power spectrum S{lo) cannot be obtained for any arbitrary switching- 
propensity function ka(t). Nevertheless, one can obtain an analytical formula for the power spectrum if the switching- 
propenstiy function is piecewise linear, as in Fig. [2]A. Fig. [2]B shows the power spectrum for a symmetric system, with 
switching-propensity functions in the forward and backward directions as shown in Fig. [^^V. It is seen that this simple, 
non-Markovian two-state model can capture the main features of the power spectrum as measured by Korobkova et 
alM. 
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III. HYBRID STOCHASTIC ALGORITHM 



The equations-of- motion for the rotor and the flagellum, Eqs. 1 and 7 of the main text, respectively, and Eqs. [T] 
and I 111 above, are propagated via a Heun scheme [l5[. 

The algorithm to determine when the next hopping, switching, or polymorphic transition will occur is essentially 
a kinetic Monte Carlo algorithm [l^. It is based on the observation that the survival probability S{t), i.e. the 
probability that no hopping, switching or polymorphic transition has happened after a time t after the last event, is 
given by 

5(t) = exp(-a(t)), (23) 
where a{t) is the cumulative total propensity function: 

a{t) = / dt'kTit'), (24) 



with k^lt) being the total propensity function as given by 

Ns 



krit) = K.^A^iit)) + k'"'"' {{Qjit)}) + fc™^rn'(^L(t) - Onit) - 0„O). (25) 

In practice, right after a hopping, switching or polymorphic transition, a random number, ^, between zero and one 
is drawn. The equations-of-motion of the rotor and the flagellum are then integrated together with the equation that 
describes the temporal evolution of a{t): 

Integrating Eq. [26]since the last event leads to an estimate for a{t) = dt'kT{t'). The next event then occurs after 
a time t since the last event when 

ait) > log(l/0. (27) 

The event type a, where a is either a hopping, switching, or polymorphic transition, is subsequently chosen with a 
probability as given by 

Pa(t) = fca(t)/fcT(t). (28) 
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